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Abstract 

We consider Standard Model Higgs boson production by gluon-gluon fusion in 
hadron collisions. We present a calculation of the next-to-next-to-leading order QCD 
corrections to the cross section in the H — > WW — *> lulu and H — > ZZ 41 decay 
channels. The calculation is implemented in the parton level Monte Carlo program 
HNNLO and allows us to apply arbitrary cuts on the final state leptons and the associated 
jet activity. We present selected numerical results for the signal cross section at the 
LHC, by using all the nominal cuts proposed for the forthcoming Higgs boson search. 



January 2008 



1 Introduction 



The search for the Higgs boson pQ and the study of its properties (mass, couplings, decay widths) 
are at the heart of the LHC physics program. In this paper we consider the production of the 
Standard Model (SM) Higgs boson by the gluon fusion mechanism. 

The gluon fusion process gg — > H, through a heavy-quark (mainly, top-quark) loop, is the main 
production mechanism of the SM Higgs boson H at hadron colliders. When combined with the 
decay channels H — > 77, H —> WW and H — ► ZZ, this production mechanism is one of the most 
important for Higgs boson searches and studies over the entire range, 100 GeV < M# < 1 TeV, of 
Higgs boson mass M H to be investigated at the LHC [21 E]. 

The dynamics of the gluon fusion mechanism in controlled by strong interactions. Detailed 
studies of the effect of QCD radiative corrections are thus necessary to obtain accurate theoretical 
predictions. 

At leading order (LO) in QCD perturbation theory, the cross section is proportional to ag, 
«s being the QCD coupling. The QCD radiative corrections to the total cross section have been 
computed at the next-to-leading order (NLO) in Refs. [U El [6] and found to enhance the cross 
section by about 80 — 100%. In recent years also the next-to-next-to-leading order (NNLO) 
corrections [3 El [H QUI dU [12] have been computed. The NNLO effect is moderate and, for a light 
Higgs, it increases the NLO cross result by about 15 — 20%. The effects of a jet veto on the total 
cross section has also been studied up to NNLO [13] . We recall that all the NNLO results have 
been obtained by using the large- M t approximation, M t being the mass of the top quark. 

The NNLO results mentioned above are certainly important, but they refer to situations where 
the experimental cuts are either ignored (as in the case of the total cross section) or taken into 
account only in simplified cases (as in the case of the jet vetoed cross section). Generally speaking, 
the impact of higher-order corrections may be strongly dependent on the details of the applied 
cuts and also the shape of the distributions is typically affected by these details. 

The first NNLO calculation that fully takes into account experimental cuts was reported in 
Ref. [14] . in the case of the decay mode H — » 77. In Ref. [T5] the calculation was extended to the 
decay mode H — > WW — *> lulu. The calculations of Refs. [T4~l[T5] were performed with the method 
described in Ref. [16], based on sector decomposition [T7] . Besides Higgs boson production, the 
above method has been applied to the NNLO QCD calculations of e + e~ — ► 2 jets [18], vector 
boson production in hadron collisions [19], and to the NNLO QED calculation of the electron 
energy spectrum in muon decay [2D] . 

In Ref. [21] we have presented an independent NNLO calculation of the Higgs production 
cross section, including the decay H — > 77. The method is completely different from that used in 
Refs. [141, [T5] . Our calculation is based on the subtraction method. 

The subtraction method [22] is probably the most popular technique to handle and cancel 
infrared singularities in QCD computations at high energy, and has lead to the formulation of 
general algorithms [231 [21] to perform NLO calculations in a relatively straightforward manner, 
once the relevant amplitudes are available. In recent years, several research groups have been 
working to develop general NNLO extensions of the subtraction method [2S1 [251 [2H [2S1 I2"9"] . 
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NNLO results, however, have been obtained only in some specific processes. The calculation of 
e + e~ — > 2 jets [301 [31] was the first to be addressed, and, more recently, the computation of 
e + e~ — > 3 jets [321 [331 El] has been completed. 

The version of the subtraction method proposed in Ref. [21] can be applied to a specific 
class of processes, namely, the production of colourless high-mass systems (lepton pairs, vector 
bosons, Higgs bosons, . . . ) in hadron collisions. As usual for calculations performed within the 
subtraction formalism, the computation can be organized into a parton level event generator. The 
latter feature is particularly useful, since the user can apply the required cuts on the final state 
and plot the corresponding distributions in the form of bin histograms. 

In Ref. [21] we have applied our method to the computation of the Higgs production cross 
section, including the decay H — ► 77. In the present paper we extend the calculation of Ref. [21] 
to the other important decay modes of the Higgs boson, namely, H —>■ WW — > lulu and H — > 
ZZ — > 4 leptons, and present predictions for the Higgs boson signal that take into account all the 
realistic experimental cuts on the final state leptons and the associated jet activity. 

The paper is organized as follows. In Sect. [2] we describe our NNLO Monte Carlo program. 
In Sect. [3] we present the results of our calculation for the decay modes H — > WW — > lulu and 
H — > ZZ —> 41. In Sect. @]we summarize our results. 



2 The HNNLO Monte Carlo program 

The numerical program HNNLO is a fortran code that implements the version of the subtraction 
method proposed in Ref. [21]. The program computes the Higgs boson production cross section 
at hadron colliders up to NNLO in QCD perturbation theory. 

The cross section up to (N)NLO can be written as 



d&(N)NLO 



n 



H 

(N)NLO 



da? n + 



LO 



H+jets CT 
a<7 (N)LO aa (N)LO 



(1) 



The first term (virtual) is the simplest to compute numerically: it contains the LO cross section 
do^o & t It = 0, Qt being the transverse momentum of the Higgs boson, suitably convoluted with a 
hard function 7i which includes the regularized one-loop (two-loop) corrections to the LO process. 
The second term (real) is the most cumbersome to evaluate. Its first contribution, da^~^, is 
the (N)LO cross section for the production of the Higgs boson in association with one (or more) 
jets. This contribution is evaluated with the version of the subtraction method of Ref. [21], as 
implemented in the MCFM [35] package. When q T — * 0, da^X^o * s divergent, and is supplemented 
with the subtraction of a suitable counterterm, da?^ LQ . The difference in the square bracket of 
Eq. ([!]) is thus finite as qr — > 0. 



In the present version of the code (version 1.1) we have implemented three decay modes for 
the Higgs boson: H — > 77 [21], H — > WW — > lulvy} and H — > ZZ 4 leptons. In the latter case 
the user can choose between H — > ZZ fi + fi~e + e~ and H — > ZZ e + e~e + e~ , which includes 



^Results for this decay channel were presented at the Les Houches Workshop "Physics at TeV Colliders" in june 
2007, and at the Radcor Conference in October 2007. 
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the appropriate interference contribution. The program can be downloaded from [SB], together 
with some accompanying notes. 



3 Results up to NNLO 

3.1 Preliminaries 

We consider Higgs boson production at the LHC (e.g. pp collisions at y/s = 14 TeV). We use 
MRST2004 parton distributions [37] , with densities and as evaluated at each corresponding order 
(i.e., we use (n + l)-loop as & t N n LO, with n — 0, 1, 2). Unless stated otherwise, renormalization 
and factorization scales are set to their default values, fiR = fip = Mh- We remind the reader 
that the calculation is done in the M t — > oo limit. As for the electroweak couplings, we use 
the scheme where the input parameters are Gp, Mz, My/ and a(M z )- In particular we take 
G F = 1.16639 x 1(T 5 GeV~ 2 , M z = 91.188 GeV, M w = 80.419 GeV and a{M z ) = 1/128.89. The 
decay matrix elements are implemented at Born level, i.e., radiative corrections are completely 
neglected^. The Higgs boson is treated in the narrow-width approximation, but in the W and Z 
decays we take into account finite width effects, by using IV — 2.06 GeV and T z = 2.49 GeV. As 
far as jets are concerned, we use the fc^-algorithm [39] with jet size D = 0.4. 

3.2 H WW -> Ivlv 

We consider the production of a Higgs boson with mass Mh = 165 GeV. The width is computed 
with the program HDECAY [40] to be Th = 0.255 GeV. With this choice of Mh the Higgs boson 
decays almost entirely into WW pairs. We consider the decay W — > Iv by assuming only one final 
state lepton combination. The corresponding inclusive cross sections are given in Table [TJ The 
NLO and NNLO i^-factors are 1.84 and 2.21, respectively, and are in good agreement with the 
inclusive i^-factors from the calculation of the total NLO and NNLO cross section [TOl [TT| fl~2]. 



a (fb) 


LO 


NLO 


NNLO 


A*f = m = M H /2 


136.37 ±0.09 


241.59 ±0.43 


268.7 ± 1.8 


= Hr = M h 


112.08 ±0.07 


206.46 ±0.33 


247.2 ± 1.3 


= = 2M H 


92.88 ±0.06 


178.43 ±0.25 


227.4 ±0.8 



Table 1: Cross sections for pp —> H + X — > WW + X — > lulu + X at the LHC when no cuts are 
applied. 

We first apply a set of preselection cuts taken from the study of Ref. [4T] . 

1. The event should contain two leptons of opposite charge having p? larger than 20 GeV and 
rapidity \y\ < 2; 

*We note that the full QCD+EW corrections to the decay modes H — > WW(ZZ) — > 4 leptons have been 
recently computed [38] . 
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2. The missing p? of the event should be larger than 20 GeV; 

3. The invariant mass of the charged leptons should be smaller than 80 GeV; 

4. The azimuthal separation of the charged leptons in the transverse plane (A</>) should be 
smaller than 135°. 

The first cut selects dilepton events originating from the decay oiWoiZ bosons. Lepton pairs 
originating from the inclusive production of a Z boson are mostly rejected with cuts 2-4. The 
corresponding cross sections are given in Table El Comparing with Table [T] we find that the 



a (fb) 


LO 


NLO 


NNLO 


Hf = Vr = M h /2 


64.03 ± 0.06 


113.57±0.28 


124.75 ± 1.28 


^F = /iR = M H 


53.10 ±0.05 


97.30 ±0.21 


116.24 ±0.81 




44.32 ± 0.04 


84.69 ±0.16 


106.48 ±0.61 



Table 2: Cross sections for pp 
cuts are applied. 



H + X — > WW + X — > lulu + X at the LHC when preselection 



efficiency is 47% both at NLO and at NNLO. The corresponding NLO and NNLO i^-factors are 
1.83 and 2.19. With respect to the inclusive case, we notice that the preselection cuts do not alter 
significantly the convergence of the perturbative expansion. 

For each event, we classify the transverse momenta of the charged leptons according to their 
minimum and maximum value, prmin and prmax- In Fig. [T]we plot the prmm and prmax distribution 
at LO, NLO and NNLO. We see that QCD corrections tend to make the distributions harder. 
This can be also appreciated from Fig. El where we compare the NNLO distributions with the 
NLO ones, normalized to the same area. 
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Figure 1: Transverse momentum spectra of the charged leptons forpp — > H + X — > WW + X — > 
lulu ± X at LO (dots), NLO (dashes) and NNLO (solid). Preselection cuts are applied. 

In Fig. [3] we plot the A(p distribution at LO, NLO and NNLO. As is well known [22], for 
the Higgs boson signal the leptons tend to be close in angle, and thus most of the events are 
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Figure 2: As in Fig. [IJ- comparison of p? spectra at NNLO (solid) with NLO normalized to the 
same area (dashes). 

concentrated at small A<f>. We notice that the steepness of the distribution increases when going 
from LO to NLO and from NLO to NNLO. As a consequence, the efficiency of a cut on this 
variable also increases with the perturbative order. 

We finally consider the following selection cuts |41j . which are designed to isolate the Higgs 
boson signal: 

1. The two charged leptons, with rapidity \y\ < 2, should fulfil prmin > 25 GeV and 
35 GeV < p rmax < 50 GeV; 

2. The missing of the event should be larger than 20 GeV; 

3. The invariant mass of the charged leptons should be smaller than 35 GeV; 

4. The azimuthal separation of the charged leptons in the transverse plane (A0) should be 
smaller than 45°; 

5. Finally, there should be no jets with p 3 ^, larger than a given value p^ to . 

These cuts further exploit: i) the shape of the PTmin and prmax distributions shown in Fig. [TJ 
ii) the strong angular correlations of the charged leptons leading to the steep A0 distribution in 
Fig. El Hi) the fact that the decay of top quarks from the tt background produces 6-jets with large 
transverse momentum. A jet veto is thus very efficient to suppress this background. 

In Table |3] we report the corresponding cross sections in the case of p£ to = 30 GeV. 

A comparison with Table [2] reveals that the cross section is strongly suppressed with respect 
to the case in which only preselection cuts are applied: the efficiency turns out to be 8% at NLO 
and 6% at NNLO. The scale dependence of the result is strongly reduced at NNLO, being of 
the order of the error from the numerical integration. The impact of higher order corrections is 
also drastically changed. The i^-factor is now 1.19 at NLO and 1.11 at NNLO. As expected, the 
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Figure 3: Normalized distribution in the variable A0 when preselection cuts are applied. 



tr(fb) 


LO 


NLO 


NNLO 


Hf = Hr = M h /2 


17.36 ±0.02 


18.11 ±0.08 


15.70 ±0.32 


= Vr = M h 


14.39 ±0.02 


17.07 ±0.06 


15.99 ±0.23 


Hf = Hr = 2M h 


12.00 ±0.02 


15.94 ±0.05 


15.68 ±0.20 



Table 3: Cross sections for pp —> H + X —> WW + X —>■ lulu + X at the LHC when selection cuts 
are applied and p™ to = 30 Ge V. 



jet veto tends to stabilize the perturbative expansion. The latter point has a simple qualitative 
explanation [T3] . 

It is well known that the effect of higher order contributions to the inclusive Higgs production 
cross section is large. The dominant part of this effect is due to soft and virtual contributions. The 
characteristic scale of the highest transverse momentum p™ ax of the accompanying jets is indeed 
^max ^ _ z ^Mh, where z = M H /s and (1 — z) measures the average distance from the partonic 
threshold. As a consequence, the effect of the jet veto is small unless p™ to is substantially smaller 
than p™ ax . Decreasing p™ to , the enhancement of the inclusive cross section due to soft-radiation 
at higher orders is reduced, and the jet veto improves the convergence of the perturbative series. 
Note, however, that when p^ to is much smaller than the characteristic scale p™ ax ~ (1 — z)Mh, 
the coefficients of the perturbative series contain logarithmically enhanced contributions that may 
invalidate the convergence of the fixed order expansion. 

In order to estimate the perturbative uncertainties affecting our calculation, in Fig. H] we report 
the LO, NLO and NNLO bands as a function of p^ eto , when all the other selection cuts are applied. 
The bands are obtained by varying [ip = hr between Mh/2 and 2M#. The results of Fig. @] deserve 
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Figure 4: Cross sections as a function of p™ to when selection cuts are applied. The bands are 
obtained by varying Hr = \if between Mh/2 and 2Mh- 

some discussion. 

At LO there are no jets accompanying the Higgs boson, and thus the cross section is inde- 
pendent on p^ to . The NLO band overlaps with the LO one for p™ to smaller than about 50 GeV. 
Without jet veto (p™ to — > oo) the i^-factor, defined with respect to the LO cross section at cen- 
tral values of the scales, ranges between 1.32 (fip = fiR = 2Mh) and 1.63 (fip = fiR = M#/2). 
Comparing with the inclusive results, we see that the selection cuts 1-4 alone already imply a 
reduction of the impact of higher order corrections. We also observe that the NLO band becomes 
very narrow as soon as p™ to decreases. 

The NNLO band overlaps with the NLO one for p^ to > 30 GeV and thus suggests a good 
convergence of the perturbative expansion in this region of p^ to . On the contrary, for p™ to <30 
GeV, the NNLO band is very narrow and does not overlap with the NLO one, suggesting that, in 
this region, the perturbative uncertainty obtained through scale variations is likely to be under- 
estimated. 

The NNLO corrections to the pp -> H + X -> WW + X -> lulu + X at the LHC were 
independently computed in Ref . [T5] . The preselection cuts we use are the same as those considered 
in Ref. [TB]. Taking into account the different normalization^], the ensuing cross sections in Table 
[2] are in good agreement with those given in Table 2 of Ref. [15] . When selection cuts are applied, 
a direct comparison is not possible, since the cuts we employ are not exactly the same. Fig. 1 
of Ref. [15| shows that, when only the jet veto is applied, the NLO and NNLO bands computed 
as in Fig. [5] overlap for p^ to < 40 GeV. Nonetheless, when all the selection cuts are applied and 
pveto _ 25 GeV, the NLO and NNLO results reported in Table 3 of Ref. [15J do not overlap. 
Although the selection cuts we use are not exactly the same, the latter result is consistent with 

Hn our calculation we strictly apply the large-M t approximation, whereas in the calculation of Ref. [TS] the 
results are normalized to the Born cross section with exact top-quark mass dependence. 
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the behaviour we observe in Fig. |H In the recent study of Ref. [43] the efficiencies obtained at 
NNLO are shown to be in good agreement with those predicted by the MC@NLO event generator 

m- 

3.3 H -> ZZ -> e + e-e + e- 

We now consider the production of a Higgs boson with mass Mh = 200 GeV. The width is 
computed with the program HDECAY [10] to be Th = 1-43 GeV. In this mass region the dominant 
decay mode is H — > ZZ —>■ 41, providing a clean four lepton signature. In the following we consider 
the decay of the Higgs boson in two identical lepton pairs. When no cuts are applied, the signal 
cross sections are reported in Table HI We find that the interference contribution is smaller than 
1% in this mass region. The ensuing inclusive cross section is thus a factor of 2 smaller than the 
cross section in the decay channel H — > ZZ — > /i + /i~e + e~Q. 



a (fb) 


LO 


NLO 


NNLO 


A*f = Hr = M H /2 


2.457 ±0.001 


4.387 ±0.006 


4.90 ±0.03 


= Hr = M h 


2.000 ±0.001 


3.738 ± 0.004 


4.52 ±0.02 


Hf = Hr = 2M h 


1.642 ±0.001 


3.227 ±0.003 


4.14 ±0.01 



Table 4: Cross sections for pp —> H + X — > ZZ + X — ► e + e e + e ± X at the LHC when no cuts 
are applied. 

The NLO i^-factor is K = 1.87 whereas at NNLO we have K = 2.26. These results are in good 
agreement with those obtained from the calculation of the total NLO and NNLO cross section 

psnniEj. 

We consider the following cuts [3]: 

1. For each event, we order the transverse momenta of the leptons from the largest (pn) to 
the smallest {pta)- They are required to fulfil the following thresholds: 

p T1 > 30 GeV p T2 > 25 GeV p T3 > 15 GeV p T4 > 7 GeV ; 

2. Leptons should be central: \y\ < 2.5; 

3. Leptons should be isolated: the total transverse energy in a cone of radius 0.2 around 
each lepton should fulfil Et < 0.05 pr] 

4. For each possible e + e~ pair, the closest (mi) and next-to-closest (m^) to Mz are found. 
Then m 1 and m 2 are required to be 81 GeV < mi < 101 GeV and 40 GeV < m 2 < 110 
GeV. 

These cuts are designed to maximize the statistical significance for an early discovery, but to keep 
the possibility for a more detailed analysis of the properties of the Higgs boson. The corresponding 
cross sections are reported in Table [51 

^In the case of H — ► ZZ — * e + e~e + e~ there is an additional diagram, obtained for example by exchanging the 
momenta of the two electrons, but there is also a symmetry factor 1/4, due to the two pairs of identical particles 

E3- 



s 



a (fb) 


LO 


NLO 


NNLO 


= Hr = M h /2 


1.541 ±0.002 


2.764 ±0.005 


3.013 ±0.023 


A*f = = M H 


1.264 ±0.001 


2.360 ±0.003 


2.805 ±0.015 


l*F = Vr = 2M H 


1.047 ±0.001 


2.044 ±0.003 


2.585 ±0.010 



Table 5: Cross sections for pp — > H + X — > ZZ + X — > e + e e + e + X at the LHC when cuts are 
applied. 

Comparing with Table HI we see that, contrary to what happens in the H — > WW — > lulu 
decay mode, the cuts are quite mild, the efficiency being 63% at NLO and 62% at NNLO. The 
NLO and NNLO if -factors are 1.87 and 2.22, respectively. Comparing with the inclusive case, we 
conclude that these cuts do not change significantly the impact of QCD radiative corrections. We 
also find that the effect of lepton isolation is mild: at NNLO it reduces the accepted cross section 
by about 4%. 

In Fig. [5] we plot the pr spectra of the final state leptons. We note that at LO, without cuts, 
the pti and pT2 are kinematically bounded by Mh/2, whereas pr3 < Mh/3 and pta < Mh/^- 
It is well known that, in the vicinity of kinematical boundaries, QCD cross sections may develop 
perturbative instabilities beyond a given order, if the behaviour of the cross section is not smooth 
at that order [IB]. This is what can be observed in the px spectra of the photons in the H — > 77 
decay mode [21]. In the present case, the effect of the cuts further reduces the kinematically 
allowed region, but the LO distributions smoothly reach their kinematical boundary, and we do 
not observe such perturbative instabilities beyond LO. 

As in Fig. [H in Fig. [5] we see that QCD corrections tend to make the distributions harder. 
This can be also appreciated from Fig. [61 where we compare the NNLO distributions with the 
NLO ones, normalized to the same area. 



4 Summary 

We have presented a calculation of the NNLO cross section for Higgs boson production at the 
LHC, in the decay modes H —>■ WW — > lulu and H — > ZZ — > 4 leptons. The calculation takes 
into account all the experimental cuts designed to isolate the Higgs boson signal [31 [JT] . In the 
case of the decay mode H — > WW — > lulu, we confirm previous findings that the effect of radiative 
corrections is strongly reduced by the selection cuts. In the case of the decay mode H — > ZZ — > 4 
leptons, we find that the proposed cuts are mild and do not change dramatically the size of QCD 
radiative corrections. 

Our calculation is implemented in the numerical program HNNL0 [3S]- The present version 
of the program includes the most relevant decay modes of the Higgs boson, namely, H — > 77, 
H — > WW — > lulu and H — > ZZ — > 4 leptons. In the latter case it is possible to choose between 
H — > ZZ — > fi + u,~e + e~ and H — > ZZ — > e + e~e + e~ , which includes the appropriate interference 
contribution. The user can apply all the required cuts on the final state leptons (photons) and 
the associated jets and plot the corresponding distributions in the form of bin histograms. These 
features should make our program a useful tool for Higgs studies at the Tevatron and the LHC. 
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Figure 5: Tranverse momentum spectra of the final state leptons for pp — > H + X — > ZZ + X — > 
e + e~e + e~ + X, ordered according to decreasing pt , at LO (dotted), NLO (dashed), NNLO (solid). 
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